Overlap Distribution of the Three-Dimensional Ising Model 
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We study the Parisi overlap probability density Ph(q) for the three-dimensional Ising ferromagnet 
by means of Monte Carlo (MC) simulations. At the critical point Pl(q) is peaked around q — 
in contrast with the double peaked magnetic probability density. We give particular attention to 
the tails of the overlap distribution at the critical point, which we control over up to 500 orders 
of magnitude by using the multi-overlap MC algorithm. Below the critical temperature interface 
tension estimates from the overlap probability density are given and their approach to the infinite 
volume limit appears to be smoother than for estimates from the magnetization. 

PACS: 05.50.+q Lattice theory and statistics (Ising, Potts, etc.), 75.40.Mg Numerical simulation 
studies, 75.10.Hk Classical spin models, 75.10.Nr Spin-glass and other random models. 
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I. INTRODUCTION 

In this paper we investigate the two replica overlap 
probability density Pl{q) for the three-dimensional (3d) 
Ising model. On a L? lattice with periodic boundary 
conditions q is defined by 
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(1) 



where and are the spins of two copies (replica) 
of the system at temperature T = 1/(3. The distribu- 
tion of the overlap q is of major importance in spin-glass 
investigations where it plays the role of an order 

parameter, often called Parisi order parameter. 

To our knowledge this quantity has never been inves- 
tigated for simple spin systems like the 3d Ising model. 
One reason is certainly that one has in that situation 
the magnetization m as an explicit order parameter at 
hand and a description of the critical properties based on 
the magnetic probability density P™ (m) is believed to be 
identical to one based on Pl(<z), in particular (q) = (m) 2 . 
However, the overlap probability density is an interesting 
object for study on its own merits and we find remark- 
able differences between the shapes of Pl(q) and P™(to). 
Therefore, we find it worthwhile to have the properties 
of Ph{q) documented for the Ising model, which is by 
orders of magnitude easier to simulate than spin glasses, 
since the dynamics is much faster and only one (instead 
of many) realization needs to be simulated. 

In the vicinity of the critical point, by finite-size scaling 
(FSS) arguments || Ph{q) can, in leading order for L 
large, be written as 



P L {q) = —P'{q') with q' = ^- 



(2) 



Here P' is a universal, L-independent function and o~l, 
is the standard deviation of q with respect to the prob- 



ability density P L (q e [-1,+1]) (or P L (q G [0, 1]) when 
appropriate) . 

A major focus of our investigation is on the tails of the 
Pl(q) distribution, which we control for L = 36 at T c 
over 500 orders of magnitude by using the multi-overlap 
MC algorithm ||. This is also of interest in view of a 
conjecture by Bramwell et al. j7|,|| that a variant of ex- 
treme order statistics describes the asymptotics of cer- 
tain probability densities for a large class of correlated 
systems. Besides the Ising model with some T(L) — > T c 
as L — > co, their class includes the 2d XY model in the 
low temperature phase, turbulent flow problems, percola- 
tion models and some self-organized critical phenomena. 
For large L the asymptotic behavior is claimed to be de- 
scribed by an L-independent curve, which for the overlap 
variable would read (q' — > oo) 



P'(q') = C exp a (q' 
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ob(<3'-<4ax) 



(3) 



Here C, a, b are constants and (/^ax = w/^Li where 
<7max is the position of the maximum of the probability 
density Pl{q) at positive q. Equation (0) is a variant of 
Gumbel's first asymptote sec r p~o|JlT| for reviews of 
extreme order statistics. 

However, Eq. (^) is in contradiction with the widely 
accepted large-deviation behavior, based on the propor- 
tionality of the entropy with the volume p2[ , 



P L (q) cx exp [-N f(q)} 



(4) 



where, for large N, f(q) does not depend on N. Our data 
support Eq. (||) . Using the multimagnetical approach |l^] 
a similar study of the tails could be performed for the 
magnetic probability density P^ l (m) 1 but this is outside 
the scope of our present paper. 

We like to point out that for the overlap distribution 
of spin glasses the status of Eq. (^) is unclear due to the 
quenched average. Our previously reported result Jl~4f] 
demonstrates that for the 3d Edwards- Anderson Ising 
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spin glass a probability distribution of the form (g) gives 
an excellent description of the tails of the Parisi overlap 
distribution. Because of the special nature of its phase 
transition, Eq. (^) may also be questioned for the 2d XY 
model, where the extreme order asymptotics (^|) (more 
precisely a variant of it) with a = tt/2 is found in the spin 
wave approximation However, the range of validity 
of this perturbative argument is unclear, at least to us. It 
may be worthwhile to employ the methods of Ref. |l3| , or 
those of the present paper, to perform a careful numeri- 
cal investigation of the the 2d XY model with respect to 
these questions. 

We have performed simulations at and below the criti- 
cal (Curie) temperature T c of the Ising model phase tran- 
sition. We approximate T c by the value of Ref. Ell 
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T'c. 



0.221654 



(5) 



and present our results for T c in Sect. |fi|. Besides address- 
ing the question of the asymptotic behavior of the overlap 
distribution, we estimate the critical exponent ratio 2/3 jv 
from the FSS behavior of the standard deviation a^. For 
the temperatures below the Curie temperature we choose 



/3 a = 0.232 and (3 2 = 0.3 



(6) 



For (3 = 0.232 multimagnetical results are available p6[ , 
which determine the probability density P™(m) of the 
magnetization over many orders of magnitude. 

Our numerical results were obtained with the spin- 
glass code of the investigations of P |l4||l7| ] by simply 
choosing all the exchange coupling constants to be equal 
to +1. A code which is specialized to the Ising model 
would be far more efficient. Therefore, we have limited 
our present simulations to small and medium sized lat- 
tices. As the results appear already quite clear, there 
seems to be no particularly strong reason to push on to- 
wards (much) larger systems. 



II. RESULTS AT THE CRITICAL 
TEMPERATURE 

Figure □ shows our overlap probability density re- 
sults PlW at the critical temperature (||). They rely 
on a statistics of 32 independent runs (with different 
pseudo random number sequences) for lattices up to size 
L = 30 and on 16 independent runs for our largest lat- 
tice, L = 36. After calculating the multi-overlap param- 
eters || the following numbers of sweeps were performed 
per repetition (i.e. independent run) : 2 19 for L — 4, 2 21 
for L = 6, 2 22 for L = 8, 2 23 for L = 12, 16, 2 24 for 
L = 24 , 2 25 for L = 30 and, again, 2 24 for L = 36 (with 
the present computer program this lattice size became 
too time consuming to scale its CPU time properly). Not 
to overload Fig. [I] error bars are only shown for selected 
values of q, whereas the lines are drawn from all data. 
The probability densities are normalized to 
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FIG. 1. Overlap probability densities at the critical point 
(3 C = 1/T C = 0.221654. 



dqP L {q) 



(7) 



and we show only the q > part because of the symmetry 
Ph{—q) = Pl(q)- We cut the range at q = 0.5, because 
for L > 8 the probability densities are almost zero for 
q > 0.5. 

Somewhat surprisingly we find the maximum of our 
Pl{q) probabilities at <7 max = 0, in contrast to the mag- 
netization where one finds a double peak at T c , see for 
instance numerical data in Ref. [ fl8f and analytical re- 
sults of Ref. |TJ| , both with periodic boundary conditions. 
In our simulation we kept a time series for the magne- 
tization, which reproduces the expected double peaked 
histograms at T c (as the accuracy of our magnetization 
histogram is lower than that of results in the literature, 
we abstain from giving a figure). Such differences are 
expected. For instance, while the two low-temperature 
magnetization values m = ±1 correspond to four overlap 
configurations, two with q = +1 and two with q = — 1, 
the inverse is not true. There are altogether 2 N overlap 
configurations with q = +1 and another 2 N with q = — 1. 

Figure || shows In [Pl(<7)] versus q. The ordinate is cut 
off at —1000, to cover a range with results from at least 
two lattice sizes. The L = 36 lattice continues to exhibit 
accurate results down to —1200, thus the data from this 
system cover 1200/ln(10) = 521 orders of magnitude. 

The collapse of the Pl (q) functions (Q) on one universal 
curve P'(q') is depicted in Fig. | The fi gure shows some 
scaling violations, which become rather small from L > 
24 onwards. The standard deviation <jl behaves with L 
according to 
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oc L- 2 ^ (l + c 2 L-" + ...) 



(8) 



Note that the ratio f}/v is defined for the magnetiza- 
tion, and by FSS theory || af oc L~^l v holds for the 
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FIG. 2. Logarithm of the overlap probability densities at 
Pc = l/T c = 0.221654. 
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sities P'(q') = o"l Pl{q) versus g' = q/cri at the critical point. 
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FIG. 3. Rescaled overlap probability densities P'(q) = 
ol Ph(q) versus q = q/cri at the critical point. 



standard deviation of the magnetization. The factor of 
two difference in the exponent of Eq. (|^) comes from di- 
mensionality. Scaling relations and estimates of the Ising 
model critical exponents are reviewed in Ref. f20|| . In 
particular, 2(3/v = d — 2 + r] holds. Our estimate of 
2(3 jv from a four-parameter fit of Eq. (ft) to our data 
for the standard deviation ctl is 2/3/ z> = 1.0293 (28) with 
Q = 0.31 the goodness of fit (see [^l[ for the definition 
of Q). Restricted to our L > 24 lattices the more stable 
two-parameter fit to the leading behavior of Eq. (||) gives 

^j- = 1.030 ± 0.005 with Q = 0.36. (9) 

The most accurate estimates of the literature |2(J clus- 
ter around r\ = 0.036 with an error of a few units in 
the last digit. Within the conventional statistical un- 



certainties this is consistent with our 2(3 /v values. The 
two-parameter fit becomes quickly inconsistent when the 
smaller lattices with L < 24 are included, with a trend 
towards smaller values of 2(3 jv. Therefore, we conjecture 
that there will be a slightly increasing trend when larger 
lattices should become available. Because its larger er- 
ror bar reflects to some extent systematic uncertainties, 
we prefer Eq. (^) over the four-parameter fit as our final 
estimate. 



TABLE I. Deviation points for A 36 ln[P' L (q')] = 0.5 of the 
L = 4 to 30 lattices from the L = 36 result. 



L 


4 


6 


8 


12 


16 


24 


30 


q 


2.27 


2.46 


2.65 


2.95 


3.25 


3.83 


4.54 


q 


0.781 


0.580 


0.472 


0.351 


0.288 


0.225 


0.212 



In Fig. ^ we show the logarithm ln[P' (<?')] of the 
rescaled overlap probability densities and we see a break- 
down of scaling for sufficiently large q' . The ordering of 
the lattices is that the rightmost curve corresponds to 
the L — 36 lattice. The smaller lattices deviate from it. 
From the left: First, the L = 4 lattice (not visible), next 
the L = 6, then L = 8, L = 12, L = 16, L = 24 and last 
L = 30. The agreement is over a larger and larger range 
in q' . However, scaled back to q, it concerns the vicinity 
of q = 0. To quantify this, we have collected in Table || 
the q' and corresponding q values at which the deviation 

A 36 HPL(Q')}^HPm(Q')}-HPL(Q')} 

becomes 1/2, a deviation too small to be visible on the 
scale of Fig. [|. The q' values are seen to increase, whereas 
the corresponding q values decrease. 

It is well known that the requirement of consistency 
of a universal probability density (Q) with the functional 
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FIG. 5. The functions extracted from Eq. (^) for 

various lattice sizes, are plotted together with a fit according 
to Eq. @. 
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FIG. 6. The functions ft(q), as well as the fit to Eq. (pr[), 
with faaiq) subtracted (error bars are indicated at selected 
values of q) . 



form (H|) determines the function /(g). Namely, to lead- 
ing order the scaling of the function (||) P' L (q r ) implies 
that 

L d / (V L-W) = g(q') (10) 

is an L-independent function. Therefore, 

f(q) cx q d »W (11) 

holds. Note that the non-critical Gaussian behavior is a 
special case, obtained for dv/2(5 = 2. 

In contrast to Eqs. (||) and (pdj), the functional form 
(|j) is expected to hold for all q, when L becomes large. 
This is easily tested by plotting 

fUq) = ~ln[P L (q)) , (12) 

as is done in Fig. ^, and seeing if /l(<?) is L-independent 
up to 0(1/N) terms as it should. Not to obscure the be- 
havior by too large symbols, the lines are plotted without 
error bars. A fit of the scaling form ( |i~lj ) to our fz§(q) 
data for q < 0.2, f(q) = 0.000049 - 0.073 q dv l 2<i with 
2/3/v = 1.030 from Eq. (^), is also included in the figure. 

We see excellent convergence towards an L-indepen- 
dent function, where the higher lying curves correspond 
to the smaller lattices (L = 4 being the one on top). 
However, the scaling behavior (|ll|) only holds in the 
vicinity of q = 0. To make this quantitatively more 
precise, we subtract the function f3e(q) from the oth- 
ers and plot the difference in Fig. ^| (at selected values 
of q we now include barely visible error bars by plot- 
ting f L (q) ± Af L (q) - / 3 6(g))- In the large volume limit 
(for both of two lattices) the difference |/li(<z) — /l 2 (<7)I 
should be bounded by a constant 




We see in Fig. ^ that the L = 24 and 30 curves fall almost 
on the L = 36 one (the zero line). Note that the figure 
is cut off at q = 0.62. The number of sweeps needed 
to propagate the system over the full admissible q-range 
scales in the multi-overlap ensemble at least proportional 
to the system size N. Aiming at a comparable statis- 
tics for all system sizes, the required computer time thus 
grows at least proportional to N 2 . Therefore and because 
of numerical problems with the floating point representa- 
tion caused by the extreme smallness of Ph{q) for q — > 1 
when L is large, we restricted the overlap simulations to 
q G [-0.7, +0.7] for the lattice sizes L = 16, 24, and 30, 
and to g 6 [-0.62, +0.62] for the L = 36 lattice. Nev- 
ertheless the smallest values of Pl(q) we sampled were 
those of the L = 36 lattice. 

On the basis of Eq. (^) the plots of Figs. || and || had to 
be expected. The conjecture (||) of Bramwell et al. [0J|] 
appears to be ruled out for the Ising model. Namely, 
when Eq. (|) (with 9l ' nax = 0) and Eq. (||) are both 
valid in some region of bq' = bL^Z", one finds to leading 
order 

bqL 2!3 ' v = d\n{L)+\n[f{q)] , (13) 

and f(q) can only be L-independent if b is not a constant, 
but depends on q and L. 



III. RESULTS BELOW THE CRITICAL 
TEMPERATURE 

Below the critical temperature we made 16 indepen- 
dent runs per lattice size with the following numbers of 
sweeps per repetition: 2 16 for L = 4, 2 17 for L = 6, 2 18 
for L = 8, 2 19 for L = 12, 2 20 for L = 16 {fi = 0.232) 
and 2 21 for L = 16 (f3 = 0.3). The overlap probability 
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FIG. 7. Overlap probability densities at P = 0.232 (left set 
of curves) and j3 = 0.3 (right set of curves). 
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FIG. 9. Logarithms of the overlap probability densities at 
(3 = 0.3. 




FIG. 8. Logarithms of the overlap probability densities at 
(3 = 0.232. 



densities Pl{<i) for f3 = 0.232 and (3 = 0.3 are shown 
together in Fig. 0. Clearly, the peaks moved away from 
zero and are now at q max — 0.3408 (L = 16, /3 = 0.232) 
and q max = 0.8237 (L = 16,0= 0.3). In Figs. | and § 
we show the logarithms of these probability densities at 
(3 = 0.232 and (3 — 0.3, respectively. The scales in these 
two figures are chosen to accommodate all the -Pl(O) 
data, but not their tails, which continue down to much 
lower values. For the largest (L = 16) difference between 
Pl (0) and the maximum of Pl (q) , we see that it increases 
from about four orders of magnitude at f3 = 0.232 to 
about 65 orders of magnitude at (3 = 0.3. 

Note that the most likely Pl(0) configurations are 
those where one replica stays around magnetization to — 
and the other around the maximum of the magnetic 
probability density at positive or negative magnetization 



TABLE II. Effective interface tension (IIS) results, F 3i l, 
from the overlap parameter at /3 = 0.232 and f3 — 0.3. At 
(3 = 0.232 results for the same quantity obtained in Ref. 
from the magnetization density are included for comparison. 



L 


= 0.232 


P = 0.232 Ref. [ L6| 


(3 = 0.3 


4 


0.00962 (27) 


0.05297 (29) 


0.23490 (33) 


6 


0.01416 (12) 


0.03403 (15) 


0.26325 (20) 


8 


0.016740 (82) 


0.02779 (13) 


0.28457 (20) 


12 


0.020281 (64) 


0.02485 (12) 


0.29751 (17) 


16 


0.022715 (34) 


0.02521 (11) 


0.29959 (11) 



to. It follows that the ratio 



Rl 



Ph (<?max) 



(14) 



where -Pl (ftnax) is the maximum of P/,(g), is related to 
the interface tension F s according to the formula intro- 
duced in CHl for the ratio PV 



Rl 



L l (0)/P£ l (qn 
CL p exp [-2L 2 F S ] + . . 



(15) 



Here C, p are constants and p — —1/2 in the one-loop 
capillary wave approximation |22| or one- loop ^-the- 
ory |T^,g3|,Q , compare the discussion in Ref. [p5| . Two- 
loop $ 4 -theory is considered in Ref. [^6|. The correction 
is large and it appears that a reliable estimate of p does 
not exist. 

To determine the interface tension one may first calcu- 
late the lattice size dependent effective interface tensions 



F, 



s,L 



lnRL ' 



(16) 



and then make an extrapolation of F s .l for L — > oo. 
Table |] collects our F s l results where the error bars 
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with respect to the last digits are given in parentheses. 
For the sake of easy comparison, we list also some F s ^l 
results of Ref. M] at f3 = 0.232, obtained by applying 
the definitions (|14| ) and ([l6]) to the probability density of 
the magnetization. It is notable that the F B l estimates 
from the overlap densities increase monotonically in the 
listed range of lattice sizes, whereas the F s> l estimates 
from the magnetization show a more complex behavior: 
Up to L = 12 they decrease, then they turn around to 
increase and the increase has been followed Jll| up to 
lattices of size L = 32. 

We pursue a similar fitting strategy as in Ref. jl6| . As 
there, it turns out that our data do not really support fits 
to more than two parameters and that including the cap- 
illary wave term with the one-loop theoretical coefficient 
p does not lead to any improvements of the goodness Q 
of the fits. In essence we are left with fits to the leading, 
likely effective, correction 



F 8 , L = F a+ a 4 



(17) 



and, due to our small lattice sizes, finite-size corrections 
are so big that the best we can do is a fit of the interface 
tensions from the largest two lattices, L = 12 and 16. 
This yields the estimates 



F s = 0.03002 (24) at f3 = 0.232 



(18) 



appears to be confined to a small q (but still large q' = 
qL 2 P/ v ) neighborhood. It may be worthwhile to check 
whether the magnetic probability distribution, for which 
comparable simulations are easier to perform, exhibits a 
similar behavior. Further, we find support in favor of 
standard large deviations @), instead of the form (|l3| ) 
derived from Gumbel's first asymptote (|J). 



Below the critical point, in Sect. Ill 



we estimate inter- 
face free energies from our overlap probability densities. 
The results are smoother than those from the probability 
density of the magnetizatio n [|16| and tend to reconcile 
discrepancies noted in Ref. [|29||. But, like at the criti- 
cal point, considerably larger lattices would be needed to 
reach high precision results. 
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and 



F s = 0.30583 (68) at f3 = 0.3 , 



(19) 



which are (under the circumstances of our limited sys- 
tem sizes) in reasonably good agreement with results of 
Hasenbusch and Pinn (HP) [^7| , for a review see Ref. [^8| . 
Again, our error bars are purely statistical and do not re- 
flect systematic errors due to our small lattice sizes. 

Our result (|l^) at (3 = 0.232 is lower than the mul- 
timagnetical estimate of Ref. |[6]]. This moves into the 
right direction and indicates that the resolution of the 
inconsistency between the multicanonical and the HP es- 
timate at f3 — 0.232, discussed in the paper by Zinn and 
Fisher J29], has its origin in the complex finite-size scal- 
ing behavior of F Sj £ estimates from the magnetization, 
which could be resolved by simulating larger systems. It 
is notable that this difficulty of the extrapolation appears 
to be limited to a small neighborhood of (3 = 0.232, as 
the multimagnetical F s estimates jl6| at (3 = 0.227 and 
(3 = 0.2439 are perfectly consistent with HP, see Fig. 1 
of Ref. M. 



IV. SUMMARY AND CONCLUSIONS 

In Sect. we have analyzed the critical behavior of 
the overlap variable q. In essence agreement with the 
standard scaling picture is found, but with some new 
insights. In particular, we exhibit in Table | that scaling 
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